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METHOD AND APPARATUS FOR DETECTING SHOCK ABSORBER DAMAGE 
BACKGROUND AND SUMMARY OF THE INVENTION 



[0001] This application claims the priority of German patent 
document 100 28 749.2, filed 10 June 2000, the disclosure of 
which is expressly incorporated by reference herein. 

[0002] The invention relates to a method and apparatus for 

detecting shock absorber damage, particularly in a motor vehicle. 

[0003] Proper maintenance of vehicles requires continuous 

availability of information concerning the wear condition of 
important vehicle components, including in particular the shock 
absorbers . 

[0004] Depending on the shock absorber construction (for 

example, single-tube gas pressure shock absorbers and twin-tube 
shock absorbers) , different damage patterns may occur, such as 
pitting on the piston rod, wear of the piston rod or leaks in the 
separating piston packing, which become noticeable by a drop of 
the shock absorbing action. The shock absorbing effect will 
decrease and the vehicle body vibrations and vehicle vibrations 
in the vertical direction will be absorbed to a lesser degree. 



[0005] In known methods for detecting shock absorber damage, 
the shock absorber is either tested in the installed condition 
on a test stand (such as the Boge Shock Tester) or is analyzed 
in the removed condition on a testing machine (such as the VDA 
Testing Machine) . These known shock absorber damage detection 
methods have the disadvantage that they provide no information 
concerning the shock absorbing characteristics during driving 
operation. Rather, time-consuming and cost-intensive checking 
is required on a special test stand; and the shock absorber may 
even have to be removed before the checking. 

[0006] One object of the invention is to provide a method and 
apparatus for detection of shock absorber damage which permit 
analysis of the shock absorber condition during driving 
operation . 

[0007] This and other objects and advantages are achieved by 
the shock absorber evaluation method and apparatus according to 
the invention, based on the idea of using the signal of an 
antilock braking system (ABS) rotational wheel speed sensor to 
detect shock absorber damage. According to the invention, a 
conclusion with respect to shock absorber condition can be drawn 
by analyzing the ABS rotational wheel speed signal in a 
particular frequency range or ranges. 
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[0008] According to a first embodiment, the rotational speed 
change An of the rim {high-pass-filtered ABS signal) is used for 
detecting shock absorber damage. 

[0009] As an alternative, the radius change Ar of the tire is 
computed from the rotational speed signal supplied by the ABS 
rotational wheel speed sensor according to the following 
equation: 

Ar = — - r 0 (1) 
n 

wherein : 

v: longitudinal vehicle velocity [m/s], if slip « 0; 

(preferably averaged from 4 rotational wheel speed signals) 
n: rotational speed of rim [wheel/s] (ABS signal, input 

quantity) 
r 0 : static tire radius [m] 

[0010] From the time varying signal for the radius change Ar 
of the tire or the rotational speed change An of the rim, its 
auto power density spectrum <3> flr or $ An is computed. The computing 
of the auto power density spectrum from the tire radius change 
signal has the advantage that it uses a quantity which is 
independent of driving speed. 
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[0011] The relationship between the road irregularity h acting 
upon the wheel (and thus upon the shock absorber) and the 
rotational speed change of the rim or the radius change of the 
tire can be represented system-theoretically by means of a 
frequency response function G fln or G Ar which also includes the 
condition of the shock absorber. 

[0012] For the excitation h {that is, the road irregularity) , 
the auto power density spectrum is given by the following 
equation : 

<& h (co) = <D h (Qo) •v M - 1 -n 0 w -ta -w (2) 
wherein: 



$ h {Q 0 ): extent of irregularity (depends on the road 
construction, such as cement base, asphalt, 
cobblestone pavement, etc.) [m 3 ] 

v: longitudinal vehicle velocity [m/s] 

w: bumpiness [-], («2) 

Q 0 W : reference wavelength [rrf w ] 

co: timing circuit frequency [wheel/s] 

[0013] (Reference is made in this respect to Mitschke, M. : 



"Dynamics of Motor Vehicles", Volume B: "Vibrations", 2nd 
Edition, Springer Publishers, Berlin, 1984) . 
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[0014] Within the scope of the invention,- it has been found 

that, in the auto power density spectrum of the rotational speed 
of the rim determined by means of the ABS rotational wheel speed 
sensor or the computed radius change of the tire, there are 
ranges which depend on the shock absorption parameter d and 
ranges which are essentially independent of the shock absorption 
parameter d. (The shock absorption parameter d characterizes the 
generated shock absorbing action with respect to the relative 
speed of the shock absorber.) In the auto power density 
spectrum, according to the invention the latter represent a 
reference frequency range which preferably has the frequencies 
co 2fi in the interval of 19 to 22 Hz or of approximately 30 to 33 
Hz. According to the invention, the ranges which are dependent 
on the shock absorption parameters represent an analysis 
frequency range in the auto power density spectrum, which 
analysis frequency range preferably has frequencies co lri in the 
interval of approximately 12 to 15 Hz. 

[0015] For the auto power density spectrum of the radius 

change of the tire 4> Ar [a ltif d) or of the rotational speed change 
of the rim d ) in the analysis frequency range or in the 

analysis interval, the following therefore applies according to 
the invention: 

<l»Ar(C0 1;i ,d) = iGArCCDud)! 2 .^^) (3&) and 
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< Wfc>i.>,d) = I G A n(o)i,i d) | 2 • $ h (a) 1>i ) 



(3b) . 



[0016] For the auto power density spectrum of the radius 

change of the tire or of the rotational speed change of the rim 
in the reference frequency range or in the reference interval, 
the following therefore applies according to the invention: 

^ArCcozj) = | G Ar (a)2,i) | 2 • 0 h (co2,i) (4a) and 

*An(002j) = I G An (0) 2>i ) f 2 . <D h (c02,i) ( 4b } ^ 

[0017J By forming a quotient of equations (3a) and (4a) as 

well as (3b) and (4b), the following equations are obtained: 



DSKW Ar = 2 L 



^(#>2,z)| 



CO\,j 
®2,i J 



(5a) and 



DSKW An =2 L 



n{co\,i>d)\ 



r \ 

C0\,i 
C0 2 ,i 



(5b) 



wherein w is preferably approximately 2 . 



[0018] According to the invention, this quotient DSKW flr or 

DSKW An corresponds to a characteristic shock absorber damage 
value, which changes when the shock absorption parameter d 
changes (that is, the shock absorber deteriorates), because, as 
a result, the frequency response function G An or G Ar is changed or 
influenced. This permits the derivation of a conclusion with 
respect to the shock absorber condition. The signal, which 
corresponds to the characteristic shock absorber damage value 
DSKW Ar or DSKW An , is preferably low-pass filtered in order to 
obtain a smoothing. 

[0019] The rotational speed signal of the rim or the computed 
radius change of the tire contains the useful signal which is of 
interest (with the recognizable resonance frequencies of the 
vehicle body and the rim or tire as well as the resonance 
frequency of the rotational speed of the rim and the belt) on the 
one hand, and a signal which is interfering for the analysis in 
the frequency range, in the form of an equal component and a 
drift in the signal. The equal component depends on the 
instantaneous driving speed, while the drift is caused by 
acceleration or deceleration of the vehicle. These interference 
effects can be eliminated, preferably by high-pass filtering of 
the rotational speed signal, particularly of the ABS rotational 
wheel speed signal, preferably with a corner frequency of 1 Hz. 



[0020] In the described embodiment according to the invention, 
computation of the characteristic shock absorber damage value 
DSKW flr or DSKW ta is a function of the parameter w (bumpiness) , 
which is preferably equal to 2. In a second embodiment according 
to the invention, a computation of the characteristic shock 
absorber damage value is suggested which is independent of the 
bumpiness, that is, of the road surface. 

[0021] Here, the basic idea is to introduce a second reference 
interval with the frequency values co 3/i , preferably in the range 
of from 30 to 33 Hz, the first reference interval co 2fi in this 
case preferably being in the range of from 19 to 22 Hz. 

[0022] The following applies with respect to the auto power 

density spectrum of the radius change of the tire or of the 
rotational speed change of the rim in a reference frequency range 
2 or in the reference interval 2: 



(6a) and 



<I>An((%i) = I G An (0)3,i) I 2 • 3>h(C03,i) 

(6b) . 



[0023] Based on the above equations (2), (3) and (4), a new 

characteristic damage value DSKW' ar or DSKW' to can now be defined. 



DSKW' Ar =2 



= 2 



= 2 



wherein : 

G x {<Oi,i, d) : frequency response function, analysis 

interval, center frequency coj, partial 
frequencies co a<i are within the analysis 
interval; 

G x (co 2(i ) : frequency response function, reference 

interval 1, center frequency co 2 , partial 
frequencies co 2 ,i are within the reference 
interval 1; 

G x {co 3/i ) : frequency response function, reference 

interval 2, center frequency co 3 , partial 
frequencies o> 3fi are within the reference 
interval 2; 

*x (oh,!, d) : auto power density spectrum radius change or 
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rotational speed, analysis interval, center 
frequency 

*x(w2,i) : auto power density spectrum radius change or 

rotational speed, reference interval 1, center 
frequency co 2 ; 

<3> x (u 3 ,i): auto power density spectrum radius change or 

rotational speed, reference interval 2, center 
frequency oo 3 ; 
angular frequency (2 • n • f) analysis interval; 

»2,ii angular frequency (2 • n • f) reference 

interval 1; 

«3,i' angular frequency (2 • n • f) reference 

interval 2; 
w : bumpiness road 

[0024] If the widths of the three intervals are now selected 
to be small, and the center frequencies a lf co 2 , and co 3 such that 
the following applies: 




• with i= l...k 



a characteristic damage value DSKW' flr or DSKW'^ is obtained which 
is independent of the bumpiness w. 
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[0025] The following is obtained for this characteristic 

damage value: 



DSKW Ar =2 



= 2 



(8a) 



DSKW' An =l 



gJ^)\ 



(8b) 



[0026] Equations (5a), (5b), (8a), (8b) show that a check has 
to take place as to whether the respective expression in the 
denominator exceeds a threshold to be defined, so that no 
division will occur by zero or by a figure close to zero. If the 
denominator expression falls below this threshold, the 
computation is not carried out. 
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[0027] Figure 11 shows the nomenclature for the frequency 

values Q Xri/ wherein x = 1, 2 or 3, in the analysis and reference 
range on the frequency axis. 

[0028] Based on equation (5a) and equation (5b), it is 

illustrated that, for the characteristic shock absorber damage 
value DSKW Ar or DSKW to , a quotient is computed from two auto power 
density values. Since the auto power density is a measurement 
of the energy content of a signal for a specific frequency 
interval, according to the invention, by means of a suitable 
mounting of filters, a characteristic damage value DSKW ArrF , 
DSKW to , F or DSKW' Ar , F , DSKW' AnrF can be defined which is similar to 
the DSKW ar/ DSKW to or DSKW' flr , DSKW T An . 

[0029] The characteristic damage value DSKW Ar , F or DSKW an , F is 

based on equation (5a) or equation (5b) . 

[0030] The definition for the auto power density spectrum is 
supplied in equation (A) . 



2T 
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[0031] This illustrates that the auto power density spectrum 
of a quantity x is equal to the squared amplitude spectrum X(jco) 
divided by twice the observation time T. 

[0032] The amplitude of a signal in a defined frequency range 
can be determined by filtering with a narrowly limited band pass, 
and subsequent value formation of the determined maxima and 
minima. The amplitude of the time signal corresponds 
approximately to the value in the amplitude spectrum for the 
corresponding frequency range. Since the auto power density 
spectrum is determined for a time period T, a filtering takes 
place by means of a first low pass (forming the average over the 
time period) . This filtering is followed by a quotient formation 
in order to obtain the DSKW ar;F or DSKW to;F (compare Figure 6), 
after further filtering by means of a second low pass (to 
determine a long-term trend) . 

[0033] Because the DSKW is a relative characteristic damage 
value, squaring of the approximated amplitude spectrum X (jco) and 
the limit value formation T - °° will not be necessary in the 
computation. Furthermore, if a constant bumpiness w (preferably 
=2) is assumed, the right-hand parenthetical expression in 
equation (5) can be eliminated. (It causes only a scaling of 
DSKW Ar or DSKW to .) The described structure for a characteristic 
shock absorber damage value DSKW ar , F or DSKW An , F is illustrated in 
Figure 6. 
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[0034] Other objects, advantages and novel features of the 

present invention will become apparent from the following 
detailed description of the invention when considered in 
conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0035] Figure 1 shows the auto power density spectrum of the 
radius change of a front tire, left; 

[0036] Figure 2 shows the auto power density spectrum of the 
radius change of a rear tire, left; 

[0037] Figure 3 shows the time variation of the characteristic 
shock absorber damage value for a front tire, left, for different 
shock absorber stages; 

[0038] Figure 4 shows the time variation of the characteristic 
shock absorber damage value for a rear tire, left, for different 
shock absorber stages; 

[0039] Figure 5 shows the position of the analysis interval, 
of the first and second reference intervals in the auto power 
density spectrum of the radius change for a rear tire, left; 
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[0040] Figure 6 is a conceptual block diagram which shows the 
computation of DSKW An , F or DSKW Ar , F ; 

[0041] Figure 7 shows the DSKW An , F course of the wheels, front 
left and rear left, for the shock absorber positions "sporty" and 
"comfortable"; 

[0042] Figure 8 is a conceptual block diagram which shows the 
computation of DSKW' Ar , F or DSKW' AnrF ; 

[0043] Figure 9 shows the analysis logic for analyzing the 

characteristic shock absorber damage values DSKW x XjX ; 

[0044] Figure 10 shows the threshold value query for detecting 
the damage; and 

[0045] Figure 11 shows the nomenclature for the frequency 

values co X/i . 

DETAILED DESCRIPTION OF THE DRAWINGS 

[0046] The auto power density spectra or characteristic shock 
absorber damage values illustrated in Figures 1 to 4 represent 
results of test drives at the BMW test site by using a BMW 740iA 
with varying shock absorbing actions (parameters d, in the 
preceding equations) . In this case, the "sporty" adjustment 
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(high shock absorption) and the "comfortable" adjustment (low 
shock absorption) were used. 

[0047 J The reciprocal values of the rotational speeds of the 
rim n from the ABS sensor signal were recorded as input data for 
shock absorber damage detection according to the invention. The 
temporal course of the radius change Ar was determined therefrom 
according to equation (1). 

[0048] As the analysis interval for the determination of the 
term <S> flr {u> lrif d) in equation (5a), the range [12-15] Hz in the auto 
power density spectrum was selected; as the reference interval 
for the determination of the term <J> Ar ( u 2,i)f the range [30-33] Hz 
was selected. The quotients determined for these frequencies 
were added according equation (5a) and were emitted after low- 
pass-filtering. 

[0049] The auto power density spectra of the radius change 3> Ar 
of the left front wheel and of the left rear wheel of the data 
used to compute the characteristic shock absorber damage value 
are illustrated in Figures 1 and 2. In the range of the analysis 
interval [12-15] Hz, the influence of the shock absorber 
condition (here "comfortable" or "sporty") upon the course of the 
auto power density spectrum is clearly recognizable. In the 
range of the reference interval [30-33] Hz, such dependence 
cannot be detected. This means that the auto power density 
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spectrum of the radius change is suitable for judging the shock 
absorber condition. 

[0050] Based on the illustrated auto power density spectra for 
the radius change, the courses of the characteristic shock 
absorber damage value for the left front wheel or the left rear 
wheel illustrated respectively in Figures 3 and 4 are obtained 
as a function of the shock absorber stage (that is, simulated 
shock absorber condition) "sporty" or "comfortable". The 
temporal courses of the characteristic shock absorber damage 
values correspond to a constant speed of 80 km/h and an identical 
route profile. A comparison of Figures 3 and 4 shows that the 
difference in the characteristic shock absorber damage value as 
a function of the selected shock absorber stage is even more 
pronounced on the rear axle. 

[0051] Various advantages are connected with the method 

according to the invention and with the system according to the 
invention. To detect shock absorber damage according to the 
invention, no additional sensor system is required, so that a 
cost-effective solution is provided by integration into an 
existing control unit. Furthermore, according to the invention, 
shock absorber condition can be analyzed during actual driving 
operation. In addition, maintenance of the shock absorbers can 
be implemented which meets the requirements. For example, shock 
absorber condition data can be provided to the respective 
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monitoring organization (such as the Technical Surveillance 
Group) during a general inspection so that the vehicle owner will 
have no additional cost for a shock absorber inspection which may 
be planned for the future. 

[0052] Using the structure illustrated in Figure 6, 

measurements were carried out at the BMW test site by means of 
a BMW 740iA. Different shock absorber actions were adjusted by 
the selection of the shock absorber position ( "sporty" 
corresponds to a high shock absorption and "comfortable" 
corresponds to a low shock absorption) . 

[0053] The pertaining DSKW An , F courses are illustrated in 

Figure 7. The top diagram shows the DSKW fln/F course for the left 
front wheel, while the bottom diagram of Figure 7 shows the 
course for the left rear wheel, once again for the respective 
shock absorber positions "sporty" and "comfortable". In partial 
areas, the algorithm is masked out, as described above. The 
illustration shows that the DSKW to , F or the DSKW flr , F is suitable for 
drawing a conclusion with respect to the shock absorber 
condition. 

[0054] The above-mentioned characteristic damage value 

DSKW' ar/F or DSKW' to>F is based on equation (8a) or equation (8b) . 
For this purpose, Figure 8 shows a simplified structure for 
determining these characteristic damage values according to the 
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invention. The signal representing the rotational speed change 
or the radius change is first filtered by means of a 
corresponding band pass filter for an analysis interval or for 
first and second reference intervals. Subsequently, a value 
formation takes place of the determined maxima and minima, 
followed by filtering by a first low pass filter. This filtering 
is then followed by a quotient formation according to equations 
8a and 8b. Finally, another filtering takes place by means of 
a second low pass filter. 

[0055] Tests using a BMW 740iA at the BMW test site were 

carried out also by means of this algorithm. Analogous to the 
above-explained approach, different shock absorber actions were 
set as a result of the two shock absorber stages "sporty" (high 
shock absorbing forces) and "comfortable" (low shock absorbing 
forces) . Qualitatively, the courses are the same as those 
illustrated in Figure 7. 

[0056] Another additional processing of the DSKW X X/X is 

illustrated in Figure 9. Here, the DSKW x ZfX can be corrected, for 
example, by way of a functional or analytical relationship or by 
an experimentally determined relationship (for example, 
characteristic diagram) with respect to driving speed. In 
addition to driving speed, corrections are also advantageous here 
by means of the throttle valve angle, torque, gear position, 
rotational engine speed and operating condition of the converter 
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clutch (in vehicles with an automatic transmission) or of the 
clutch switch (in vehicles with a manual transmission) . By 
selecting the corner frequency of the "low pass 3" filter, the 
characteristic shock absorber damage value DSKW x XrX can be 
averaged over a long term (for the useful life of the shock 
absorber) . 

[0057] Shock absorber damage can be detected by querying the 
threshold value. As soon as the uncorrected or corrected DSKW X X x 
is in the damage range (Figure 10), a conclusion is drawn that 
a shock absorber is defective and this information is processed. 
On the one hand, the information concerning a shock absorber 
defect can immediately be transmitted to the shop (defect memory 
with read-out in the shop or active shop notification) ; and on 
the other hand, it can be made accessible to the customer in the 
form of a display concept. 

[0058] The foregoing disclosure has been set forth merely to 
illustrate the invention and is not intended to be limiting. 
Since modifications of the disclosed embodiments incorporating 
the spirit and substance of the invention may occur to persons 
skilled in the art, the invention should be construed to include 
everything within the scope of the appended claims and 
equivalents thereof. 
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